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ABSTRACT 

The statistical property of the weak lensing fields is studied quantitatively 
using the ray-tracing simulations. Motivated by the empirical lognormal model 
that characterizes the probability distribution function of the three-dimensional 
mass distribution excellently, we critically investigate the validity of lognormal 
model in the weak lensing statistics. Assuming that the convergence field, k, 
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is approximately described by the lognormal distribution, we present analytic 
formulae of convergence for the one-point probability distribution function (PDF) 
and the Minkowski functionals. The validity of lognormal models is checked in 
detail by comparing those predictions with ray-tracing simulations in various cold 
dark matter models. We find that the one-point lognormal PDF can describe the 
non-Gaussian tails of convergence fields accurately up to v ~ 10, where v is 
the level threshold given by v = ^/(/t 2 ) 1 / 2 , although the systematic deviation 
from lognormal prediction becomes manifest at higher source redshift and larger 
smoothing scales. The lognormal formulae for Minkowski functionals also fit 
to the simulation results when the source redshift is low, z s — 1. Accuracy 
of the lognormal-fit remains good even at the small angular scales 2' 9 ^ 4', 
where the perturbation formulae by Edgeworth expansion break down. On the 
other hand, the lognormal model enables us to predict the higher-order moments, 
i.e., skewness, S$ jK and kurtosis, S^ K , and we thus discuss the consistency by 
comparing the predictions with the simulation results. Since these statistics 
are very sensitive to the high(low)-convergence tails, the lognormal prediction 
does not provide a quantitative successful fit. We therefore conclude that the 
empirical lognormal model of the convergence field is safely applicable as a useful 
cosmological tool, as long as we are concerned with the non-Gaussianity of v £ 5 
for \ow-z s samples. 

Subject headings: cosmology: theory — gravitational lensing — large-scale struc- 
ture of universe — methods: numerical 



1. Introduction 

Cosmic shear, coherent distortions in galaxy images caused by the gravitational field of 
the intervening large-scale structure, has now been recognized as a powerful cosmological tool 
(see Mellier 1999 and Bartelmann & Schneider 2001 for reviews). Since the signal of cosmic 
shear reflects the gravitational potential of the total mass distribution, the cosmic shear can 
be a direct probe of the dark matter distribution as well as a way to measure the cosmological 
parameters. Recently, measurements of the cosmic shear have been independently performed 
by several groups (Van Waerbeke et al. 2000a; Wittman et al. 2000; Bacon, Refregier & 
Ellis 2000; Kaiser, Wilson & Luppino 2000; Maoli et al. 2001; Rhodes et al. 2001; van 
Waerbeke et al. 2001a). The estimated shear variance from the different observational data 
set quantitatively reconciles with each other, which is in good agreement with theoretical 
predictions based on the cluster normalized cold dark matter (CDM) model of structure 
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formation. Ongoing and future wide-field surveys using high-resolution CCD camera will 
promise to dramatically improve the signal-to-noise ratio of the weak lensing signals. Hence, 
a detail understanding of the weak lensing statistics will be necessary to extract all the 
cosmological information present in the data. 

Of course, simplest but useful statistical measure of the weak lensing field is popular 
two-point statistics (Blandford et al. 1991; Miralda-Escude 1991; Kaiser 1992), however, a 
more clear discriminator is needed to break the degeneracy in constraining the cosmological 
parameters. In this respect, the non-Gaussian features of weak lensing field can be very 
sensitive to the density of universe and the higher-order statistics such as the skewness of 
local convergence field are expected to constrain the density parameter Qq and even to the 
cosmological constant A (Bernardeau, Van Waerbeke & Mellier 1997; Jain & Seljak 1997; 
Van Waerbeke et al. 2001b). On the other hand, topological analysis using the Minkowski 
functionals of the convergence field has been proposed by Sato, Takada, Jing & Futamase 
(2001) in order to extract the non-Gaussian features (see also Matsubara & Jain 2001). 
Unfortunately, their methodology heavily relies on the validity of the perturbation theory 
of structure formation and cannot be applied to the weak lensing field at small angular 
scales 9 4', where the underlying three-dimensional mass distribution is in the highly 
nonlinear regime. While the non-Gaussian signature of weak lensing fields becomes easily 
measurable on small scales, the non-Gaussian signal primarily reflects the nonlinear growth 
of the density fields, which strongly depends on the dark matter clustering properties. It is 
therefore desirable to investigate the statistical properties of three-dimensional mass density 
and explore the relation between mass distribution and weak lensing field. 

Very recently, Kayo, Taruya & Suto (2001) performed a detailed analysis of the one- and 
two-point statistics of the three-dimensional mass distribution using the high-resolution N- 
body simulations. They found that the lognormal models of one- and two-point probability 
distribution functions (PDF) can provide an excellent approximation to the nonlinear mass 
distribution with Gaussian initial condition, irrespective of the shape of the initial power 
spectrum. The lognormal distribution has been long known as an empirical prescription 
characterizing the dark matter distribution and/or the observed galaxies (e.g., Hamilton 
1985; Bouchet et al. 1993; Coles & Jones 1991; Kofman et al. 1994), however, there exists 
no rigorous explanation for its physical origin. Nevertheless, the lognormal model is now 
widely utilized in the astrophysical context such as analytical modeling of dark halo biasing 
and Ly-a forest (e.g., Mo & White 1996; Taruya & Suto 2000; Bi & Davidsen 1997). In this 
regard, the result of Kayo et al. (2001) is interesting and can be useful in quantifying the 
dark matter distribution. Furthermore, their results indicate that the lognormal distribution 
can also describe the weak lensing field, since the weak lensing effect primarily reflects the 
clustering property of dark matter distributions. Then, crucial but fundamental questions 
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arise as follows. How does the lognormal property emerge in the weak lensing field ? Can 
the lognormal model provide an approximate and reliable prescription for the weak lensing 
statistics ? 

In this paper, we quantitatively investigate those issues using the ray-tracing simula- 
tions. We first consider how the statistical feature of three-dimensional mass distribution 
is related to that of the local convergence, i.e, the two-dimensional projected density field. 
Assuming that the convergence field is well-approximated by the lognormal distribution, we 
derive analytic formulae for the one-point PDF and the Minkowski functionals. We then 
perform quantitative comparisons between the lognormal models and the results of ray- 
tracing simulations. Furthermore, we discuss the consistency of the lognormal model with 
the higher-order statistics on which the previous works have mainly focused. 

The plan of paper is as follows. In section 2, we briefly describe the basic definitions 
of weak lenses. Then, we discuss the statistical relations between mass density field and 
convergence field in section 3. The analytic expressions for the one-point PDF and the 
Minkowski functionals are presented. In section 4, the detailed comparison between the 
lognormal predictions and the ray-tracing simulations is described. The consistency between 
lognormal model and previous study is discussed in section 5. Finally, section 6 is devoted 
to conclusions and discussion. 



2. Preliminaries 

Inhomogeneous mass distribution of the large-scale structure deflects the light ray tra- 
jectory emitted from an angular direction 6 in the source plane by a small angle 56 in 
the image plane. The differences between the deflection angles of light rays emitted from a 
galaxy thus induce a distortion of the galaxy image characterized by the following symmetric 
matrix: 

where di represents a derivative with respect to 0j, and <fi is the gravitational potential of the 
three-dimensional mass density field. The variable Xs means the quantity x & t the source 
redshift, and the quantities x an d r(x) respectively denote the comoving parts of the radial 
and the angular-diameter distance: 

Xiti = / 77T7T ! H{z) = #oaM(1 + zf + (1 - - A )(l + zf + A , (2) 
Jo H \ z ) 
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and 

r sm(VK X )/VK (K > 0) 

r(x) = \ X (K = 0) (3) 

[ sinh( v / ^x)/v /Z ^ (K<0) 
with the quantity K being the spatial curvature of the universe, K = (H /c) 2 (Qo + Ao — 1). 
The matrix $j j is usually decomposed into the trace part, i.e, the convergence field n e and 
the tidal shear components ~f defined by 

«fl = ($1,1 + $2,2) , (4) 

le = -*2,2) -i^i,2- (5) 

The shear field 7 e can be directly estimated from the observed ellipticity of the galaxy 
images. The convergence field n e can be then reconstructed from the shear map from the 
relations (4) and (5) (e.g., Kaiser & Squires 1993). Since we are interested in the statistical 
properties of the weak lensing fields, one can safely employ the Born approximation, where 
the quantity didj<p(x) are computed along the unperturbed photon trajectory (Blandford et 
al. 1991; Miralda-Escude 1991; Kaiser 1992). From equations (1) and (4), the convergence 
field along a line-of-sight is simply expressed as a weighted projection of the mass density 
fluctuation field 5 (e.g, Mellier 1999; Bartelmann & Schneider 2001 for reviews): 

rxs 

Ke (0)= dxw(x,Xs)8[r(x)0,x], (6) 
Jo 

where the weight function w(x, Xs) is 

Note that, even in the weak lensing limit, the density fluctuation 5 is not small but can 
become much larger than unity. 

In statistical analysis of weak lensing fields, the smoothing filter is practically used in 
order to reduce the noise due to the intrinsic ellipticity of galaxies. Throughout this paper, 
we adopt the top-hat smoothing function. The variance of the local convergence can then 
be expressed as 

<* 2 > = £ d X Kx, X,)] 2 / |^ P^(k±) W* u (k ± r(x)9 m ) , (8) 

where P mass (/c) is the three-dimensional power spectrum of the mass distribution and the 
function Wth{x) represents the Fourier transform of the top-hat smoothing kernel: 

# TH (x) = 2^M. (9) 

Ob 
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Note that equation (8) is derived by employing the Limber's equation (Kaiser 1992). Simi- 
larly, variance of the gradient field, ((V/t) 2 ), can be expressed as 

<(V*0 2 > = £ S d X Hx,Xs)] 2 r\ X ) W^(k L ri X )6 m ) . (10) 

Equations (8) and (10) are valid on relevant angular scales we are interested in, where the 
small-angle approximation holds with good accuracy (Hu 2001). 



3. Lognormal model prescription 

3.1. One-point PDF 

The expression (6) gives a simple interpretation that the statistical feature of k is closely 
related to that of the underlying density field 5. To show this more explicitly, we consider 
the one-point PDF, which can be constructed from a full set of the moments. Let us write 
down the PDF of the mass density field 5: 

/ +oc dr 
-— ^ e -™(SM+Mx)/* 2 dSj (n) 

where a denotes the root-mean-square (RMS) of S, a — (5 2 ) 1 / 2 and (fs(x) means cumulant 
generating function: 

Similarly, the one-point PDF of the local convergence is written as 

P{k)dK = e -«(A/<*">^)^(,)/<A»> ^ (13) 

with the cumulant generating function for k, (pk(x). Here we have introduced the normalized 
convergence field, k = /c/|/c m in| so as to satisfy the range of the definition being — 1 < k < 
+oo. The quantity K m i n denotes the minimum value of the convergence field. Note that the 
actual value of K m [ n in the universe should, in principle, depend on the nature of dark matter 
between the source galaxies and observer (Metcalf & Silk 1999; Seljak & Holtz 1999). If the 
dark matter is composed of a compact object such as the primordial black holes, K, min mig ht 
correspond to the convergence evaluated along an empty beam (see eq.[31]). On the other 
hand, when the dark matter is composed of a smooth microscopic component as is suggested 
by the standard CDM scenarios, K min becomes larger than the empty-beam value for relevant 
smoothing scales (see also Jain, Seljak & White 1999). As will be discussed later in section 
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4.2, K min is an important quantity in controlling the lognormal model and is sensitive to the 
cosmological model. 

Even though we have the simple relation (6) between 5 and k, the relation between 
the cumulant generating functions of ips and is practically intractable without recourse 
to some assumptions or approximations. Under the hierarchical ansatz for the higher order 
moments of 5, i.e., S ri: s = const., Munshi & Jain (2000) and Valageas (2000) showed that 
the relation between ipg and ip^ is greatly reduced and is expressed in a compact form: 

<p k (x)= ^TT^^h r~r^ x ' ( 14 ) 



h(0) V l«min| {k 2 } 

where the quantity I K is defined by 

W) = J W?e(*±r(x)fl). (15) 

This is valid as long as the small-angle approximation holds. As also proposed by those 
authors, since the factor w I K /(\ K, m i n \(k 2 )) in equation (14) has typical values of order unity, 
one might use a simple approximation 

(pa(x) ~ (fs(x). (16) 

It is easy to check that this approximation holds for the limit of z s — > 0, where the both 
quantities w J K /( | /t m in | (k 2 )) and fdx(k 2 ) / I K (9) approach unity. 

Once equation (16) is given, PDF of the convergence field, P(k), can be directly calcu- 
lated from the one-point PDF of the three-dimensional mass density field P(S), irrespective 
of the projection effect (see eq.[6]): 



K 




dn 


i r a 

\ ^min | 


I ^min | / 


I ^min | 



P(k) dn= P[5^ a -> \ ' , (17) 



Now, let us recall the empirical fact that the PDF of local density P(S) is approximately 
described by the lognormal distribution: 



P " {S) dS = ^WTW) exp | 2in(i + <H) ) T+7' (18) 

Substituting (18) into (17), one finally obtains 



[ln(l + /t/|/t min |) + a 2 n /2] 2 | d_K_ 

2na hl '" r [ 2cr m J k + |/c n 



fln(«) dK = exp <| ■ i^Hi^V^jru^ ^ _^ _ (1Q) 



- 8- 



where the quantity <7 m is denned by 



2 _ 



= In 1 + 




) 



(20) 



2 



mm 



Before discussing the application and validity of equation (19), two important remarks 
should be mentioned. First, recall that the lognormal distribution violates the assumption of 
hierarchical ansatz used in the derivation of equation (14). The resultant expression (19) is 
thus inconsistent, however, the violation of the hierarchical assumption is fortunately weak 
on nonlinear scales of interest here (e.g., Kayo et al. 2001). We therefore expect that, as 
long as we restrict ourselves to the range of applicability, the expression (19) can provide 
a reasonable approximation and capture an important aspect of the non-Gaussianity in the 
weak lensing field. Secondly, notice that the result (19) heavily relies on the validity of the 
approximation (16). One might suspect that the approximation breaks down for a case with 
the high source redshift such as z s ^ 1, where the lensing projection becomes more important, 
and this will be discussed below. Of course, one can directly evaluate the convergence 
PDF P(k) from equations (13) and (14) assuming the lognormal PDF of 5, although this 
treatment is not useful in practice. Instead, in this paper we a priori assumes the simple 
analytic prediction (19), whereby we can further derive the useful analytic formulae for the 
Minkowski functional and obtain a more physical interpretation of the results. We will 
then carefully check our model by comparing the predictions with the numerical simulation 
results. 



Expression (17) implies that the projection effect is primarily unimportant and is elim- 
inated by rescaling the quantities 5 — > /t/|/t m m|- That is, the statistical feature of k directly 
reflects that of three-dimensional density 5. If this is the case, one can further develop a 
prediction of the weak lensing field based on the lognormal ansatz. To investigate this is- 
sue, the other informative statistics such as the higher-order correlation and the isodensity 
statistics should be examined. 

Among these statistics, it is known that Minkowski functionals can be useful and give 
a morphological description to the contour map of the weak lensing field (e.g, Schmalzing & 
Buchert 1997). In a two-dimensional case, Minkowski functionals are characterized by the 
three statistical quantities: the area fraction, v , circumference per unit length, v\ and Euler 
characteristics per unit area, t>2. The third functional is equivalent to the famous genus 
statistics often used in the cosmological context (Gott, Melott & Dickinson 1986). These 
quantities are evaluated for each isocontour as a function of level threshold, v = k/ {k 2 )- 



3.2. 



Minkowski functionals 
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The explicit expression of Minkowski functional vq(v) is given by 

v (is) = (0 (« - v (k 2 ) 1 ^)) , (21) 

where 6 denotes the Heaviside step function. The Minkowski functionals Vi(v) and v 2 (u) 
provide an additional information on the convergence field, since the definitions include 
higher derivative terms such as <9;K and didjK. According to Matsubara (2000), Minkowski 
functionals v\ and v 2 in the two-dimension map can be related to the level-crossing and the 
genus statistics defined in the two-dimensional surface, N 1 and G 2 (see also Schmalzing & 
Buchert 1997): 

Vl (u) = |N 1 (i/) = |(fe(«- (k 2 ) 1 / 2 !/) |^|>, (22) 
V2 (u) = G 2 (u) = (6 d (k- (k 2 ) 1/2 u) 5 D {d 1 K) \8 2 k\ d 1 d 1 K). (23) 



Assuming that the convergence field is approximately described by the lognormal distri- 
bution, the analytic expressions for Minkowski functionals can be derived in a straightforward 
manner. Since the area fraction represents the cumulative probability above the threshold 
v, substituting (19) into the definition (21) yields 



r +o ° i r 

Vo,\n(u) = dn Pin(ft) = -erfc < 

with the function y(y) given by 



V2 



(24) 



y{y) ^ £k + ln(l + ,(. 2 )VV|. min |)^ (25) 

2 (7i n 

As for the Minkowski functionals V\ and v 2 , the analytic expressions for the lognormal 
distribution can be obtained from the local transformation of Gaussian formulae, Ni(v) and 
G 2 {v). The details of the derivation are described by Taruya & Yamamoto (2001) (see also 
Matsubara & Yokoyama 1996). The final results become 

„ llnM = J_ * i^ 1 ' 2 e -*™ (26) 

Note that in the limit (/t 2 ) 1 ^ 2 /|/t m in| <C 1, the above predictions recover the Gaussian formu- 
lae: 

^oW = ^ erfc (£)> (28) 
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1/2 



Ul,o(f) 



8x/2 v <* 2 > ; 

1 <(Vk;) 2 > 




(29) 



^2, G (^) 



-^ 2 /2 



(30) 



2(2tt) 3 / 2 (k 2 ) 



Apart from the PDF P(k) and the area fraction vq, there exists no clear reason that 
the extension of lognormal model to the Minkowski functional v\ and vi still provides a 
good approximation in the real universe. At present, the predictions (26) and (27) are 
just regarded as an extrapolation from the one-point PDF P(k) and should be checked by 
numerical simulations. 



In this section, we quantitatively examine the validity of lognormal model by comparing 
its predictions with the simulation results. A brief description of ray-tracing simulation is 
presented in section 4.1. In section 4.2, the model parameters in the lognormal formulae 
are checked in details using the simulation data. Some important numerical effects are also 
discussed. Then, section 4.3 describes our main results, i.e, the comparisons of one-point 
PDF and Minkowski functionals between the lognormal predictions and the simulation data. 



In order to investigate the lognormal property of convergence fields, we use a series 
of ray-tracing simulations in three cold dark matter models (SCDM, LCDM, OCDM for 
Standard, Lambda, Open CDM models, respectively). The cosmological parameters used 
here are summarized in Table 1. 

To perform a ray-tracing simulation, a light-cone data set is first generated by particle- 
mesh (PM) N-body code. The PM code uses 256 2 x 512 particles and is performed in a 
periodic rectangular box of size (L, L, 2L) with the force mesh 256 2 x 512. The initial condi- 
tions are generated according to the transfer functions of Bond & Efstathiou (1984). Then 
the light-cone of the particles is extracted from each simulation during the run (Hamana, 
Colombi & Suto 2001). The final set of light-cone data are created so as to cover a field of 
view of 5 x 5 square degrees and the box sizes of each output are chosen so as to match to 
the convergence of the light ray bundle. As a result, the angular resolution of ray-tracing 
simulation, which is basically limited by the spatial resolution of PM N-body code, becomes 
almost constant, 9 ICS 1.5' from the observer at z = (Hamana & Mellier 2001). 



4. Comparison with ray-tracing simulations 



4.1. Ray-tracing simulations 



- 11 - 



Once the light-cone data is obtained, ray-tracing simulations are next performed using 
the multiple lens-plane algorithm (e.g., Schneider, Ehlers & Falco 1992; Jain, Seljak & White 
2000; Hamana, Martel & Futamase 2000). In our calculation, the interval between lens planes 
are fixed to comoving length, 80/i _1 Mpc. The 512 2 light rays are then traced backward from 
the observer's point to the source plane. The initial ray directions are set on 512 2 grids 
and the two-dimensional deflection potential are calculated by solving ray-bundle equations 
keeping the same grids in each lens plane. We obtained 40 realizations by randomly shifting 
the simulation boxes. 

After constructing the weak lensing map in the image plane at z — 0, the smoothed 
convergence field are finally computed on 512 2 grids by convolving the top-hat smoothing 
kernel. All the statistical quantities such as one-point PDF and Minkowski functionals are 
evaluated from the data. In a subsequent analysis, we use the convergence data set fixing 
the source redshifts to z s — 1 and 2. 



Since the lognormal predictions presented in previous section heavily rely on the three 
parameters, K m i n , (k 2 ) and ((V/t) 2 ), we first check them in some details using the ray-tracing 
simulations. 

Figure 1 shows the minimum value of the convergence, K min against the smoothing 
angle. The error-bars indicate the 1-a errors around each mean, where the mean value K min 
is obtained from 40 realization data. As stated in section 3.1, the convergence k evaluated 
along an empty beam is theoretically considered as a possibility of K, min and expressed from 
equation (6) as 



which implies that the light ray propagates through empty space with 5 — — 1 everywhere 
along the line of sight. In figure 1, thick lines represent the theoretical prediction (31) for 
each model. Clearly, these predictions give much smaller values than the simulation results in 
both z s = 1.0 and 2.0 cases. In addition to the systematic cosmological model dependences, 
the scale-dependence of K m \ n also appears. 

These results mean that none of the light rays becomes completely empty beam (see also 
Jain, Seljak & White 1999). Nevertheless, one cannot exclude the possibility that some or 
even the majority of lines of sights become empty at scales smaller than the mean separation 
angle between particles of N-body simulations. In this sense, the minimum value /t min could 
be affected by the finite sampling from a limiting survey size of the convergence map. Based 



4.2. On the lognormal model parameters 




(31) 
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on this consideration, we propose the following intuitive way to explain the behaviors of 
/t min shown in figure 1. Assuming that the original PDF of infinitesimal light rays obeys the 
lognormal model with the minimum value K empty , one can roughly estimate K min from the 
condition that the expectation number of independent sampling area for /t ompty < k < « m i n 
becomes unity in the field-of-view 5° x 5°. The details of estimation is described in Appendix 
A. For a given angular scale 9 TH and cosmological model, the quantity « m i n is evaluated by 
solving the equation (A2) with the prior PDF (Al). The thin lines in figure 1 depict the 
estimation based on this prescription, which reasonably agrees with simulations. This result 
is interesting in the sense that we could analytically predict « m i n and the resultant minimum 
value is sensitive to the cosmological parameters, especially Q . The practical possibility will 
be again discussed in the final section. 

Next, turn to focus on the variances of the convergence and its gradient fields. In figure 
2, the measured amplitudes, (k 2 ) 1 ^ 2 and ((V/t) 2 ) 1 / 2 , are plotted in linear- and log-scale and 
are compared with linear (short- dashed) and nonlinear (solid) predictions. For the RMS 
of the convergence (ft 2 ) 1 / 2 , the nonlinear prediction based on the Peacock & Dodds (1996) 
formula faithfully reproduces the simulation result over the smoothing angles, 9 VCS ~ 1.5'. 
On the other hand, in the case of ((V/t) 2 ), nonlinear predictions systematically deviate from 
simulations. The discrepancy remains even at larger smoothing angles, 9 ~ 10' and causes 
30 ~ 40% error. 

While the theoretical predictions based on the expressions (8) and (10) were computed 
assuming that the mass power spectrum P uiass (k) are continuous and has infinite resolution, 
the simulation data is practically affected by the finite resolution. In our case, the limitation 
of PM force mesh could be attributed to the finite resolution or the cutoff of Fourier modes, 
which becomes influential even on larger smoothing angles, 9 £ 9 TCS . In fact, compared 
to (k 2 ), dominant contribution to the quantity ((V/t) 2 ) comes from the short- wavelength 
modes of three-dimensional density fluctuations, sensitively depending on the choice of the 
smoothing filter. To show the significance of this effect, nonlinear prediction is modified 
according to the finite resolution of PM code (see Appendix B). Shortly, the Fourier- integrals 
in the expressions of (8) and (10) are discretized so as to mimic the PM N-body resolution 
as follows (Eqs.[Bl][B2] in Appendix B): 



POO 



iVn 



dk ± A 2 (k ± )W 2 n (k ± r9) J2 Ak ^ A2 ( k ±,n)W 2 u (k^ n r9), 

71=1 

( (k ± /2n) P^ ss (k±) for (/t 2 ) , 

A 2 (A:x) = 

( (kl/2n) P mass (A; ± ) for ((Vk) 2 ). 

where n-th Fourier mode k± jTl is given by k^ n = Ak± x n and the interval Ak± is set to 
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Ak±_ = 2n/L hox (z). 

The long-dashed lines in figure 2 show the results including the finite- resolution effect. In 
each panel, the number of Fourier modes, iV max is chosen as 90 (see discussion in Appendix 
B). It is apparent that the prediction of ((V/t) 2 ) 1 / 2 systematically reduces its power all 
over the scales and it almost reconciles with the simulation result, although the result of 
(ft 2 ) 1 / 2 remains unchanged. We have also examined various cases by changing the maximum 
number to 64 ^ N max ,$,128, but obtained qualitatively similar behavior. Note that the 
incorrect prediction of ((V/t) 2 ) leads to a systematic error in predicting the amplitude of the 
Minkowski functionals v\ and V2 from their definitions (26) and (27), while it does not alter 
the non-Gaussian shapes of those functionals with respect to the level threshold. 

Keeping up the above remarks in mind, in what follows, to make a comparison with 
the lognormal models transparent, we use the parameters (k 2 ), ((Vk) 2 ) and /t min directly 
estimated from the simulations when plotting the lognormal predictions (19), (24), (26) and 
(27). 

4.3. Results 

4.3.1. One-point PDF 

As a quick view of the validity of lognormal model, we first deal with the one-point PDF, 
P{k). Figures 3 and 4 show the one-point PDFs of the local convergence in various CDM 
models with the smoothing angles, 9 — 2', 4' and 8'. Here, the PDF data is constructed by 
binning the data with Ak = 0.01. The source redshifts are fixed to z s = 1.0 in figure 3 and 
z s = 2.0 in figure 4. Clearly, the functional form of the one-point PDF becomes broader as 
increasing the source redshift z s . On small angular scales, the PDF significantly deviates 
from the Gaussian PDF denoted by dashed lines. Although local convergences on larger 
smoothing scales tend to approach the Gaussian form, they still exhibit a non-Gaussian tail. 

In figures 3 and 4, solid lines represent the lognormal predictions (19), whose parameters 
« m i n and (k 2 ) are directly estimated from simulations. The agreement between the lognormal 
model and the simulation results is generally good. In particular, at lower source redshift z s = 
1.0, the lognormal PDF accurately describes the non-Gaussian tails in the high- convergence 
region up to k ^ 10 (/t 2 ) 1 / 2 . In the case of higher source redshift z s = 2.0, the discrepancy 
becomes evident at larger smoothing scale. The lognormal PDF over-predicts in the high- 
convergence region, and under-predicts in the low-convergence region. This discrepancy 
might be ascribed to the projection effect of gravitational lensing (see eq.[6]), since the 
lognormal PDF is obtained based on the approximation (17), which cannot be validated 
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as increasing the source redshift z s . Nevertheless, even in that better agreement 

between lognormal model and simulations was found at small angular scales. These features 
are also seen in figure 5, where the differences between simulated PDF and lognormal PDF 
normalized by the simulated PDF, [P S i m (/c) — -Pi n ( K )]/-P S im( K ) are plotted as a function of 
level threshold v — k/ (k 2 ) 1//2 , in the case of LCDM model. 

Accurate lognormal-fit in the low-z s case is amazing and is regarded as a considerable 
success. As a closer look at the non-Gaussian tails, however, lognormal PDF slightly under- 
predicts the simulations at low-density region, k < (see Fig.3 and upper-panel of Fig. 5). 
Furthermore, at 9 — 2', the simulation results generally tend to deviate from lognormal PDF 
in the highly non-Gaussian tails n ^ 10 (ft 2 ) 1 / 2 , although the prediction still remains consis- 
tent within the l-a error. In particular, in the case of the OCDM model, the discrepancy 
becomes apparent even at k ^ 8 (/t 2 ) 1 / 2 . The systematic deviation seen in non-Gaussian tails 
might be ascribed to the presence of virialized objects. Notice that the angular scale 9 = 2' 
corresponds to the effective smoothing scale R ~ 0.5/i _1 Mpc at mean redshift z ~ 0.5. This 
indicates that high-K value is attained by the light-ray propagating through the high-density 
region in the very massive halos with M ^ 1O 14 /i~ 1 M . Indeed, highly non-Gaussian tails 
of convergence PDF sensitively depends on the detailed structure of non-linear objects, as 
pointed out by Kruse & Schneider (2000). They construct an analytic model of one-point 
PDF based on the universal profile of dark matter halos and the Press-Schechter theory for 
halo abundance. Since the virialized halos induce the highly non-Gaussian tails of the local 
convergence, their treatment would be helpful to describe the non-Gaussian tails of PDF. 
In contrast, due to the lack of physical bases, no reliable prediction is expected from the 
empirical lognormal model. 

Another reason for discrepancy might be the choice of the parameter, K min . Strictly 
speaking, minimum value of k seen in the simulated PDF represents a rare event for the 
whole data sets, which are not rigorously equivalent to the averaged value of 40 realization 
data shown in figure 1. We have also examined the lognormal-fit adopting the minimum value 
of the PDF data for « m i n . The thin-lines in figure 5 show the results in LCDM case adopting 
the actual minimum value k of each PDF data. The results seem somehow improved at 
small angle 9 = 2'~ 4' in high but we rather recognize the fact that lognormal-fit 

to the high-density region is sensitive to the choice of K m [ n . 

Except for these details, lognormal model of one-point PDF remains a fairly accurate 
model of the convergence field, at least, up to k ~ 5(/t 2 ) 1 / 2 and is indeed applicable even at 
small angular scales such as 2' £ 9 £ 4', irrespective of the assumption, (16) or (17). 



4-3.2. Minkowski junctionals 



Having recognized the successful lognormal-fit to the one-point PDF, we next investigate 
the lognormal model of Minkowski functionals. For this purpose, we restrict our analysis to 
the source redshift z s = 1.0. Here, the Minkowski functionals for a simulated convergence 
map are calculated by the method developed by Winitzki & Kosowsky (1997). 

Figure 6 shows the results in various CDM models fixing the smoothing angle 6 = 
2.0'. The Minkowski functionals are plotted against the level threshold v = k/(k 2 ) 1 ^ 2 with 
the interval, Au = 0.5, so that data in each bin is approximately regarded as statistically 
independent. The solid lines depict the lognormal predictions (24), (26) and (27). 

Similar to the one-point PDF, a significant non-Gaussian signature is detected from 
the asymmetric shape of the Minkowski functionals, especially from the Euler characteristic 
v 2 . In marked contrast with the Gaussian predictions, the lognormal predictions remarkably 
reproduce the simulation results, not only the shape dependences but also the amplitudes. 
The agreement between lognormal prediction and simulations still remains accurate over the 
rather wider range, —4 < v < 4, where the discrepancy seen in the one-point PDF of the 
OCDM model is not observed. Since the prediction has no adjustable parameter and only 
uses the information from an output data, this agreement is successful. 

Figure 7 depicts the results with various smoothing angles fixing the cosmology to 
LCDM model. For illustrative purpose, the amplitudes of v\ and v 2 at the smoothing angle 
— 4' and 8' are artificially changed by multiplying the factors as indicated in each panel, in 
order to make the non-Gaussianity manifest. The Minkowski functionals tend to approach 
the Gaussian prediction as increasing the smoothing angle, consistent with the behaviors of 
one-point PDF. The results in other cosmological models are also similar and the agreement 
between lognormal models and simulations is satisfactory. 

To manifest the accuracy of the lognormal-fit in contrast to the other existing analyt- 
ical models, let us now consider the perturbation predictions. Employing the Edgeworth 
expansion, the perturbative expressions for Minkowski functional are derived analytically 
(Matsubara 2000; Sato et al. 2001): 

voiy) * v , G (u) + (K 2 )^-^=e-^S^H 2 (u), (32) 




(33) 



o(2) 

H 2 (v) + -*-} • (34) 



-16- 



Notice that in the case of the three-dimensional density field, the above expansion is valid 
up to the RMS of 5, a £ 0.3 (Matsubara & Yokoyama 1996; Matsubara & Suto 1996), which 
can be translated into the condition, (k 2 ) 1 ^ 2 £ 0.3 |/t m i n |- Here, the function H n (y) denotes 
n-th order Hermite polynomial, H n (x) = (— l) n e x2 / 2 (d/dx) n e~ x2 / 2 , and the quantities 
represent the skewness parameters defined as 

o(0) (or,) 

5 d) _ 3 (" 2 -v 2 "> m 

bs - 4 (k 2 )((Vk) 2 )> [6b) 

(2) ((V«-V«)V 2 «) ^ 

Equations (32)- (34) imply that, in the weakly nonlinear regime, the non-Gaussian features 
on the Minkowski functionals can be completely described by the above parameters, which 
can be evaluated by the second-order perturbation theory of structure formation (e.g., 
Bernardeau et al. 1997). It is worth noting that in the case of the lognormal model, all 
the skewness parameters, are equal to 3 (Hikage, Taruya & Suto 2001). 

Figure 8 plots the perturbation results in LCDM model. The source redshift is fixed to 
z s = 1.0. The solid lines represent the results in which the skewness parameters are evaluated 
using the second-order perturbation theory (perturb 1), while the dashed lines depict the 
results using those estimated from simulations (perturb 2). In both cases, the variances (k 2 ) 
and ((Vk) 2 ) in the expressions (32)-(34) are estimated from simulations. Also, in figure 9, 
the comparison between various model predictions is summarized, introducing the fractional 
error, defined by 

(sim) / \ (model) / \ 



v- 



with the quantity v\ f l ™l x being the maximum value of Vi(u) among the mean values of the 
simulation for each Minkowski functional: f^mlx = ^mix = 0.040 and i^mix = 0.0042 
at angular scale 9 = 2', for instance. As a reference, the 1-a error of simulation results is 
plotted as error-bars around zero mean in each panel. 

At the large smoothing angle 9 = 8', both of the perturbation results tend to reconcile 
with each other and fit to the simulation results well within the 1-a error. As the smoothing 
angle decreases, however, the perturbation results cease to fit the simulations, because the 
RMS of local convergence reach (k 2 ) 1 / 2 >, 0.4 |/t min | and the Edgeworth expansion breaks 
down. Furthermore, the second-order perturbations under-predict the skewness parameters, 
compared with those estimated from simulations, which lead to the different predictions 
(compare perturb 1 with perturb 2 in left- and middle-panels in Figs.8, 9). In particular, 
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the circumference v\ shows a peculiar behavior, v\ < 0, which is not allowed by definition. 
Notice that even in these cases, the fractional error Err[t>j(z/)] in the lognormal prediction 
still remains smaller, although the systematic deviations in every model are not so large. 

From these discussions, we conclude that the empirical lognormal models can provide a 
good approximation for the Minkowski functionals, compared to the current existing models. 

5. Consistency with higher-order moments 

Since the skewness of the local convergence has been proposed as a simple statistical 
estimator of non-Gaussian signature to determine the cosmological parameters, various au- 
thors have investigated the usefulness of this quantity using ray-tracing simulations. These 
analyses have revealed that the skewness at small-angle 9 ^ 10' exhibits the significant influ- 
ence of nonlinear clustering and a more reliable theoretical model beyond the perturbation 
theory is needed. According to these facts, non-perturbative predictions based on the "hyper- 
extended perturbation theory" (Scoccimarro & Frieman 1999) or non-linear fitting formula 
of bi-spectrum(Scoccimarro & Couchman 2001) are examined (Hui, 1998; Van Waerbeke et 
al. 2001b). 

In general, the one-point PDF as well as the Minkowski functionals characterizes a 
family of higher-order statistics. Hence, as a consistency check of the lognormal prediction, it 
seems natural to analyze the higher-order moments of local convergence. From an empirical 
lognormal PDF (19), the skewness and the kurtosis of the local convergence defined by 
S 3jK = (k 3 )/(k 2 ) 2 and S^ K = ((/t 4 ) — 3(k 2 ) 2 )/ (k 2 ) 3 are respectively given by 

** = irq( 3+ irrV)' (39) 

| r^min | \ | "Trim | / 

^^{ 16+15 FT +6 (F^) 2 + (Fn 

I ""mm | | | "ram | \ | "riim \ / \\ "Trim | 

Figure 10 shows the direct measurement of Ss )K and S^ K fixing the source redshift to 
z s = 1.0. The error-bars indicates 1-a error estimated from the 40 realization data. The 
lognormal predictions (39) and (40) are depicted as solid lines. 

In practice, the convergence field k in numerical simulation does not extend the entire 
range between K m [ n and +oo, but is limited as /t rmn < k < /t max , due to the finite sampling 
effect from a limited size of simulation data (see Kayo et al. 2001 in the case of three- 
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dimensional density field S). Thus, the n-th order moments of k given by 

(K n ) = dn P hl (K) n n (41) 

may provide a better description in evaluating the skewness S$ iK and the kurtosis S^ K . 

The dashed lines in figure 10 show the lognormal prediction based on the equation 
(41), where the cutoff parameter K max is adopted as the averaged maximum value from each 
realization data. In contrast to the accurate fit seen in the one-point PDF, the lognormal- 
fit of skewness and kurtosis seems very poor. The prediction without cutoff K max tends 
to over-predict as increasing smoothing angle, which shows opposite behavior compared to 
the simulations. As for the lognormal model based on (41), while the predictions reduce 
their amplitudes and broadly agree with simulations in the case of S^ K , they still exhibit 
some systematic discrepancies in S^ K . We have also examined the lognormal-fit adopting the 
minimum and maximum values of PDF data itself, but the result is not improved drastically. 

At first glance, these discrepancies seem to contradict with the results in one-point 
PDF, however, a closer look at one-point PDF reveals that the lognormal PDFs slightly 
under-predict the simulation results at the low-density region. This tiny discrepancy may 
be ascribed to the overestimation of skewness. In other words, the skewness as well as 
kurtosis is very sensitive to the rare events, i.e, high- and low-convergence parts of the non- 
Gaussian tails. This sensitivity is clearly shown in the OCDM model. The simulated PDF at 
smoothing angle = 2' exhibits a highly non-Gaussian tails and it overshoots the lognormal 
prediction (see upper-right panel in Fig.3). The resultant skewness and kurtosis yield values 
larger than those of lognormal prediction. On the other hand, in SCDM and LCDM models, 
lognormal PDFs accurately fit to the highly non-Gaussian tails and thereby the predictions 
of skewness and kurtosis becomes relatively consistent with simulations, at least on the small 
scales, 2'< 9 < 4'. 

Therefore, the lognormal model of convergence does not provide an accurate predic- 
tion for the statistics sensitive to the rare events. This disagreement simply reflects the 
fact that the empirical lognormal model does not correctly describe projected structure of 
dark matter halos. On the other hand, a sophisticated non-perturbative model based on the 
hyper-extended perturbation theory or fitting formula of bi-spectrum is constructed so as to 
reproduce the N-body results of higher-order moments, which can provide an accurate pre- 
diction for the convergence skewness (Hui, 1998; Van Waerbeke et al. 2001b). Hence, for the 
prediction of higher-order moments, the non-perturbative model is more useful and reliable 
than the lognormal model. In contrast, the lognormal model fairly describes non-Gaussianity 
around the peak of the PDFs over the broad range, (k 2 ) 1 / 2 , which cannot be described 

by such a non-perturbative model. Thus, at least as a complementary approach, the log- 



-19- 



normal model is useful beyond the perturbative theory, and applicable in characterizing the 
one-point PDF and the Minkowski functionals. 

6. Conclusion and Discussion 

In the present paper, we have quantitatively investigated the extent to which the log- 
normal model fairly describes the statistics of weak lensing fields on linear and nonlinear 
scales using the ray-tracing simulations. The validity of lognormal model has been checked 
in details by comparing the lognormal predictions of the one-point PDF and the Minkowski 
functionals of convergence field with their simulation results. 

The convergence field seen in the one-point PDF and the Minkowski functionals displays 
the non-Gaussian feature and significantly deviates from the Gaussian predictions. We have 
shown that the analytic formulae for lognormal models are useful and accurately describe 
the simulation results on both small and large smoothing angular scales in the case of the 
\ow-z s data, while the perturbative prediction by Edgeworth expansion fails to reproduce 
the simulation results on small scales because of the nonlinearity of the underlying three- 
dimensional density field. The detailed comparison revealed that the lognormal model does 
not provide an accurate prediction for the statistics sensitive to the rare events such as 
the skewness and kurtosis of the convergence. We therefore conclude that, as long as we 
are concerned with the appropriate range of the convergence, (k 2 ) 1 / 2 , the lognormal 

model empirically but quantitatively gives a useful approximation characterizing the non- 
Gaussianity features on the convergence field. 

The results obtained here will lead to an important improvement for the estimation of 
cosmological parameters using the Minkowski functionals (Sato, Takada, Jing & Futamase 
2001). Although the original methodology has been proposed with the use of the Edgeworth 
formulae, (32)-(34), a more reliable estimation of cosmological parameters will be possible 
using the lognormal models. In the light of this, a reliable theoretical prediction for the model 
parameters of the lognormal formulae K min , (k 2 ) and ((V/t) 2 ) should be further explored. As 
shown in figures 1 and 2, the simulation results reveal that K m i n is very sensitive to the 
density parameter Q , while the (k 2 ) and ((Vk) 2 ) depend on Q and a 8 - Such predictions 
will offer a new opportunity to determine the cosmological parameters, independently of the 
method using the higher-order moment, S 3:K . In particular, if we do not focus on the overall 
normalization of the Minkowski functionals v i and t>2 controlled by the (k 2 ) and ((Vk) 2 ), the 
non-Gaussian features on those functionals are primarily controlled by /t m i n , and, for example, 
the intuitive way to predict K m { n discussed in Appendix A could be used in constraining Qq. 
Of course, for proper comparison with observation, there are other systematic effects we have 
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to take into account. One is the redshift distribution of source galaxies. Another important 
effect is the contamination by the intrinsic source ellipticity. This not only reduces the 
signal-to-noise ratio, but also systematically affects the non-Gaussianity of the weak lensing 
signal (Jain, Seljak & White 2000). These effects should be correctly incorporated into the 
lognormal predictions. These issues are now in progress, and will be presented elsewhere. 

The statistics of cosmic shear directly reflect the statistical feature of mass distribution 
and using this fact, one might even discriminate the nature of dark matter. Furthermore, the 
weak lensing statistics have a potential to reveal the nonlinear and stochastic properties of 
galaxy biasing. In any cases, we expect that the lognormal property of the weak lensing field 
is helpful and plays an important role in extracting the various cosmological information. 

We thank J.Sato for providing us the code to calculate the Minkowski functionals from 
the ray tracing simulation data. We also thank Y.Suto for careful reading of the manuscript 
and critical comments, CHikage for invaluable discussion, T.Buchert and M.Bartelmann for 
useful comments. I.K is supported by Takenaka-Ikueikai Fellowship. T.H. and M.T acknowl- 
edge supports from Japan Society for Promotion of Science (JSPS) Research Fellowships. 
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A. Influence of finite sampling on estimation of K min 

Minimum value of the convergence, « m i n is theoretically expected to be equal to /t ompty 
(eq.[31]). In practice, however, the minimum value estimated from the simulation data can 
be systematically larger than those obtained from the empty beam, due to the finite sampling 
from the limiting survey area. 

To see the influence of the finite sampling effect, let us roughly estimate the minimum 
value, K m in. Assuming that the prior one-point PDF, P (prior) (/t) is approximately described 
by the lognormal PDF, in which the minimum value is characterized by /t cmpty instead of the 
actual value 



"mm- 



P<^(«) = -=L- exp (_Nl + ^/l^ T l) + ^/2] 2 l dn 

V27TO-ln I 2(7 ln J K +Km pty | 

with the quantity af n being ln(l + (/t 2 )/|/t C m P ty| 2 )- Then the probability that the observed 
minimum value /t min systematically deviates from /t empty is given by 



l K min| Trfft 

dKP^\K) = -^, (A2) 

|lompty| field 



where # TH is the top-hat smoothing angle, #fi c id is field-of-view angle. The right hand side 
of equation (A2) represents the lower limit of probability determined from the expectation 
number of independent sampling area. The minimum value « m i n is obtained by solving 
equation (A2). 

Based on the equation (A2), the resultant values of K m i n are summarized in figure 1 
(thin-lines) . Here, the field-of-view angle #fi e id is fixed to 5° and the parameters in the prior 
lognormal PDF, (k 2 ) and /t ompty are computed according to the theoretical predictions (8) 
and (31), respectively. 



B. Effect of finite resolution and variances (k 2 ) and ((V/t) 2 ) 

Statistical analysis based on the N-body simulation with PM algorithm should be care- 
fully treated if we deal with the statistics on small scales. In our ray-tracing simulation, 
the box size of each simulation, L hox (z), are determined so as to satisfy the resolution angle 
6 VCS ~ 1.5', i.e, L hox (z) m r(x(z)) rcs . The mesh of the PM algorithm is fixed to 252 2 x 512 
in each simulation box. Among these parameters, the number of mesh severely restricts the 
Fourier modes of mass fluctuations, which could affect the evaluation of (k 2 ) and ((Vk) 2 ), 
depending on the choice of smoothing filter. 
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In order to investigate the significance of finite mesh-size, the theoretical predictions 
(8) and (10) are modified according to the PM N-body code. Since the influence of finite 
mesh-size mainly affects the high-frequency part of the fluctuations, the Fourier-integrals in 
the expressions (8) and (10) are modified to 



Jo 



dk ± A 2 (k ± )W 2 u (k ± r6) =s> A ^ & 2 (k±,n)W* H (k ± , n r6), (Bl) 

^ n=l 

where 

r (k ± /2n) P mass (A; ± ) for (k 2 ), 
A 2 (k L ) = I (B2) 

I (fci/27r) P mass (A; ± ) for ((V/t) 2 ). 

The n-th Fourier mode k± j7l is given by k± j7l = Ak± x n and the interval A/c_|_ is set to 
Ak ± = 2n/L hox (z). 

In the above modification, the number of Fourier modes, iV max might be crucial in 
evaluating the quantity sensitive to the high frequency mode, which is related to the number 
of mesh, N mcsh = 256. Recall that the Nyquist frequency restricts the high-frequency mode to 
^Nyq = (Ak±/2)N mesh , which implies iV max = iV mes h/2. Further, the number of independent 
Fourier mode is reduced by the factor 1/2 in evaluating the power spectrum. Thus, the 
high-frequency cutoff is roughly given by fc cut ~ \[2 x Ak±(N mcsh /A), which yields 

iV max = ^ 90. (B3) 

The long-dashed lines in figure 2 represent the results taking into account the finite 
mesh-size. Here, the nonlinear mass power spectrum P mass (/c) by Peacock & Dodds (1996) is 
used in evaluating the discretized Fourier- integral (Bl). We also examined the various cases 
by changing the maximum number, 64 < iV max < 128, but obtained qualitatively similar 
behavior: the prediction ((V/t) 2 ) systematically reduces its power all over the scales, while 
the amplitude of (k 2 ) almost remains unchanged. Further, similar modification to the z- 
integral has been made so as to match the number of multiple lens-plane, however, this does 
not affect the final results. 
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Table 1. Cosmological parameters used in 
N-body simulations. 



Model 
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h 


08 


SCDM 
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0.6 


LCDM 


0.3 
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0.9 


OCDM 


0.3 


0.0 


0.7 


0.85 
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Fig. 1. — Minimum values of the convergence field, ft min , evaluated from the simulations. 
Thick lines indicate the predictions assuming the empty beam (see eq.[31]). The thin lines 
represent the estimation taking account of the effect of finite sampling (see Appendix A). The 
vertical dotted lines denote the resolution limit of ray-tracing simulation, which is primarily 
determined by the force resolution of PM N-body code (see text in Sec. 4.1): z s = 1.0 (upper 
panel); z s = 2.0 (lower panel). 
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Fig. 2. — RMS of convergence field (upper-panel) and its gradient field (lower-panel) as a 
function of smoothing angle, plotted in linear- and log-scales. The error-bars indicate the 
simulation results, where error-bars are estimated from the 40 realization of ray-tracing data 
by randomly shifting each simulation boxes. The solid and short-dashed lines represent the 
theoretical prediction based on the expressions (8) and (10). In evaluating these equations, 
nonlinear mass power spectra by Peacock & Dodds (1996) are adopted in solid lines, while the 
linear power spectra are used in short-dashed lines. The long-dashed lines also represent the 
nonlinear prediction, but taking account of the finite resolution of PM N-body simulations 
(see Appendix B for details): SCDM model (Left) ; LCDM model (Middle) ; OCDM model 
(Right) of z s = 1.0 cases(Upper-panet) and z s = 2.0 cases(Lower-panet) . 
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Fig. 3. — One-point PDFs of the convergence field fixing the source redshift to z s — 1.0 with 
smoothing angle 8 = 2', 4' and &(top to bottom). Solid lines show the lognormal prediction 
(eq.[19]), where the parameters k and (k 2 ) are directly estimated from simulations. For 
comparison, the Gaussian PDF with the same variance (k 2 ) are plotted as dashed lines: 
SCDM model (left) ; LCDM model (middle) ; OCDM model (right). 
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Fig. 4. — Same as Figure 3, but we fix the source redshift z s = 2.0. 
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Fig. 5. — Differences of one-point PDFs normalized by the simulated PDF, [P S im( K ) — 
-Pin(ft)]/-Psim(ft) as a function of level threshold v = ^/(k 2 ) 1 ^ 2 in LCDM model. 
Upper(Lower)-panel shows the results fixing the source redshift to z s = 1.0 (2.0). The 
solid, short-dashed and long-dashed lines represent the cases with smoothing angle 9 = 2.0', 
4.0' and 8.0', respectively. In plotting the ratios, the lognormal prediction adopting the 
averaged minimum value K min (see Fig.l) is used in thick lines, while the thin lines represent 
the results adopting the minimum value of each PDF data. 
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Fig. 6. — Minkowski functionals as a function of level threshold v — k/ (k 2 ) at the smoothing 
angle 9 = 2.0'. The source redshift is fixed to z s = 1.0. Solid lines show the lognormal 
predictions based on the formulae (24), (26) and (27), where all the parameters are estimated 
from simulations. The dashed lines are the Gaussian predictions obtained from (28), (29) 
and (30): SCDM model (left); LCDM model (middle); OCDM model (right). 
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Fig. 7. — Same as Fig.6, but we fix the cosmological model to LCDM model. Here, for 
illustrative purpose, the amplitude of v\ and V2 at the smoothing angle 9 = 4, 8 [arcmin] are 
enhanced in order to clarify the differences as indicated in each panel: 9 = 2'' {left); 9 = 4' 
(middle); 9 = 8' (right). 
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Fig. 8. — Edgeworth expansion of Minkowski functional compared with the simulations. In 
each panel, the results of z s — 1 case in LCDM model are shown. Solid lines (perturb 1) 
indicate the perturbation results, where the skewness parameters are calculated via pertur- 
bation theory, while the skewness parameters of the dashed lines (perturb 2) are estimated 
from simulation data directly. In plotting both cases, the variances (k 2 ) and ((V/t) 2 ) are 
fitted to the numerical simulations. Note that the amplitude of v\ and i>2 at the smoothing 
angle = 4, 8 [arcmin] are enhanced in order to clarify the differences as indicated in each 
panel. 
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Fig. 9. — Fractional errors of various model predictions in Minkowski functionals, 
Err[i>j(i/)] = [f - sim) (z/) — v { ™ AA) (v)\j v { *™l^ in the case of LCDM model with z s = 1. Here, 
v i*mlt denotes the maximum value of Vi(v) estimated from simulation. In each panel, the 
fractional errors for the lognormal and Gaussian prediction are plotted as thick-solid and 
dashed lines, respectively. The dotted- and dot-dashed lines represent the perturbative pre- 
dictions based on the Edgeworth expansion (dotted: perturb 1; dot-dashed: perturb 2). In 
each panel, the error-bars around zero mean indicate the realization error of ray-tracing 
simulations: 9 = 2' (left); 9 = 4' (middle); 9 = 8' (right). 
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Fig. 10. — Skewness S 3tK (upper) and kurtosis S^ K (lower) as a function of smoothing angle. 
The crosses, hexagons and squares with error-bars represent the simulation results in SCDM, 
LCDM and OCDM, respectively. The solid lines show the lognormal prediction based on 
the expressions, (39) and (40). The dashed lines also indicate the lognormal prediction, but 
here we take into account the limited range of convergence data, i.e., K min < k < K max . 



